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Abstract 

We discuss localization of quasiparticles in an extended NS structure in the 
situation when the reflection from the NS interface is mostly of Andreev 
kind. The localization of quasiparticle states arises due to trajectory retracing 
caused by Andreev reflection. This effect is semiclassical in the sense that in 
the classical limit fi — * the states become fully localized, while quantum 
diffraction effects destroy the localization behavior. We derive the criteria 
for this localization and show that it takes place only for sufficiently smooth 
disorder with large correlation length, whereas short-range disorder destroys 
the effect. Localization of quasiparticle states gives rise to resonance peaks in 
the local tunneling density of states. 
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We consider excitations in a thin normal metal film on the surface of a superconduc- 
tor. We assume that electrons with energy well below superconducting gap A are reflected 
from the NS interface by Andreev mechanism, which is the case if the NS interface is 
clean and smooth (see Fig. |l|). If Andreev reflection occures exactly backwards |IJ, the 
electron orbit will be closed. Such an electron will bounce between metal-vacuum and 
metal-superconductor interfaces, after each two reflections returning exactly to its original 
position. In such a model, all trajectories appear to be localized. This simple fact can be 
interpreted as integrability of classical Andreev billiards @. 

Below we discuss quantum localization in this system. We add details to the oversim- 
plified picture discussed above and correct it. It turns out that in the classical Andreev 
dynamics localization (and integrability) is not a universal behavior. We discuss complica- 
tions arising in the quantum problem, such as the role of smooth disorder versus short-range 
disorder. We compare the situation with that in a perfectly uniform metal film, where the 
localization is absent. However, in presence of smooth disorder with large correlation length 
localization can take place. We derive the criterion for localization of this type and discuss 
robustness of this phenomenon. 

Recently, the effect of Andreev reflection was discussed in the context of the problem of 
quantum chaos in NS cavities |3|J3]]. It was found that in the presence of Andreev scattering, 
chaotic dynamics gives rise to an energy gap centered at the Fermi level. This gap can serve 
as a benchmark of chaotic dynamics. In contrast with the cavity problem, we study an 
infinite system. In this case, a natural equivalent of ergodic (chaotic) dynamics in cavities 
will be derealization of states over the entire system. Accordingly, in the infinite system, 
an equivalent of regular dynamics in cavities is localization of states. 

We begin the discussion by reviewing properties of the states in a uniform film of constant 
thickness. In this system electronic states are plane waves and thus are not localized. To 
see that, one can solve Bogoliubov-de Gennes equations || for this problem and derive the 
spectrum of electrons M: 
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Here n is the quantum number (n = 0, 1,2, . . .), L is the thickness of the film, p x is the 
momentum of the particle parallel to NS interface, p F is Fermi momentum, v F is Fermi 
velocity (here we set H = 1). The spectrum (|I|) has dispersion (v(e) = de/dp x ^ 0), and 
therefore the states are not localized. However, the dispersion is much weaker than for free 
electrons and thus the Andreev states are much closer to becoming localized. 

One can qualitatively understand this dispersion as follows. When Andreev reflection 
occurs, an electron is converted into a hole, and its energy (measured from e F ) changes sign. 
Therefore, the momenta of electron and hole, p e and ph, are related by 

2 2 

^- = t F + e ; = t F - e (2) 

2m ' 2m K ' 

Here is Fermi energy, and m is effective mass. From (|2]) it is clear that for the energy 
not right at the Fermi level, i. e. for e ^ 0, the momenta p e and ph are different. On the 
other hand, the component of momentum parallel to the interface is conserved at Andreev 
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reflection. Consequently, the reflection angle must change. For e <C ep the change in the 
angle is small, and the Andreev reflection law reads: 



sm6 r V €p) 

(6i and 9 r are the angles of incidence and reflection measured from the normal to the in- 
terface). Thus, after two subsequent reflections the particle will not return to the starting 
point. Instead, it will be displaced along the interface (see Fig. |l].a). The drift velocity de- 
rived from this argument is the same as the one obtained above from dispersion in (|l|). Note 
that the same effect is responsible for the suppression of spectral flow in the superconducting 
vortex core 0. 

Evidently, even though there is no localization in a uniform system, by making the 
interface rough one can reach the situation when classical Andreev trajectories will become 
localized, with no drift along the interface. Thus our next step is to introduce disorder in 
this model. 

We consider here a model of disorder in which the thickness of normal film is slowly 
varying. More precisely, we assume that the normal metal-vacuum interface has some 
roughness, whereas the NS interface is flat. The main effect of interface curvature is that 
it acts like focusing mirror, counteracting the dispersion. To overcome the dispersion, the 
curvature must exceed certain threshold which will be estimated below. 

Note that the localization effect we consider is based on the semiclassical picture of a 
particle (nearly) retracing its trajectory after being Andreev-reflected. Therefore, the scat- 
tering by the surface roughness must also be semiclassical in order to preserve the trajectory. 
The point is that quantum effects in scattering, i. e. diffraction, can destroy the localization. 
Indeed, due to diffraction, quantum scattering is stochastic and thus it violates reversibility 
of individual trajectories. (Because two scattering events on the same disorder configura- 
tion may not lead to identical results.) Thus, we consider only sufficiently smooth surface 
fluctuations and formulate below a quantitative condition on the degree of smoothness. 

Suppose that the mean thickness of the normal metal film is L, the variation of the 
thickness is A, and the spatial scale on which the thickness varies is r c (see Fig. §). The 
criterion for localization can be expressed in terms of these parameters. 

To begin with, let us ignore diffraction, and consider a purely classical motion. Without 
any loss of generality, we can limit the discussion to the problem in a two-dimensional space. 
We will use coordinate system in which the NS interface is the line y = 0, and the metal- 
vacuum interface is given by y = L(x). It is very instructive to consider a metal-vacuum 
surface of constant curvature, i. e. of a spherical shape. 

First, suppose that the NS interface is in the equator plane. Consider a trajectory 
which hits the NS interface exactly at the center of the sphere (see Fig. |3|). Note that any 
such trajectory retraces itself even if Andreev scattering does not occur exactly backwards. 
Thus, for this trajectory the drift is indeed eliminated by the curvature. One can also show 
geometrically that if the center of hemisphere is within the normal region, i. e., L > R, 
there will be no net drift (see below). Thus, the situation when the center of curvature lies 
exactly on the NS interface, is critical. 

The conventional approach to dynamical problems of this kind involves mapping of 
Poincare section |7j in the phase space. For the D = 2 problem the phase space is four- 
dimensional, but the energy conservation reduces the dimension to three. Hence, Poincare 
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section is two-dimensional. To construct it, we consider the points on the metal-vacuum 
interface hit by electron (and disregard the points due to holes). These points are character- 
ized by their ^-coordinates Xj. Also, we characterize the momentum of electron by the angle 
9i between the momentum direction and the normal to the surface at the collision point. 
Thus, each collision is gives a pair (x{, 9i), and the trajectory of the particle is represented as 
a sequence of points in the (x, 9) plane. The Poincare section for a particular metal-vacuum 
surface shape is shown in Fig. [|. This section exhibits typical Kolmogorov-Arnold-Mozer 
features: stable periodic islands representing finite motion, i. e., non-escaping trajectories, 
and the regions around the islands representing escaping trajectories. The islands correspond 
to localized states, whereas the outer regions correspond to delocalized states. 

To find the region of stability of localized orbits, one should consider the stability of a 
self-retracing trajectory, like the one in Fig. |3|. Suppose that the particle starts from metal- 
vacuum interface at a point x and its momentum direction is characterized by the angle 9 
defined above. Suppose also that <C 1. Using the reflection law ([|) and simple geometrical 
considerations, one can write down the linearized equations for the phase space coordinates 
x' and 9', describing the state of quasiparticle after after reflection from the NS interface 
and returning to the metal- vacuum interface: 

(£) = am(S) , (4) 

where a — e/ep *C 1, and 

M(a)=( a f-i l\ 1 , „ t, L \ ■ (5) 



After one Andreev reflection the sign of e changes, because electron turns into hole. Since 
in the construction of Poincare section we are interested only in electrons and not in holes, 
the matrix describing the motion in the (x, 9) plane, is M(a)M(—a). In the first order in a 
its eigenvalues are 




(6) 

To the same accuracy in a, it is more natural to write 



Ai, 2 = exp ^±2zct^l - L/Rj (7) 

to assure the phase volume conservation. The eigenvalues Ai and A2 are complex for L < R 
and real otherwise. It means that the self-retracing trajectory is stable for L > R. 

To derive this relation in a more intuitive way, note that for each orbit the maximal 
value of x corresponds to 9 = 0, i. e., to normal reflection. Attentive reader will notice 
that here we have the same situation as in the above example with the hemisphere. When 
x is maximal, which corresponds to the trajectory turning point, the trajectory is almost 
self-retracing, like that in Fig. [3[ The center of curvature at the point x lies on the NS 
interface. We conclude that localized states exist if the center of curvature is above the NS 
interface. 
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The previous discussion shows that the localization criterion is L > R, where R is the 
curvature radius. By estimating R ~ r^/A, one arrives at 

r c < r* ~ VTK . (8) 

We call this condition classical criterion of localization. 

The criterion (|8]) is not the only constraint on r c . In this system, localized and delocalized 
classical trajectories can coexist at the same energy. (Indeed, for any value of e there are 
trajectories going along straight lines parallel to the NS interface.) Hence, any perturbation 
which mixes these two types of states can destroy localization. In particular, due to finite 
size of surface fluctuations which focus electrons, quantum diffraction takes place. Efect of 
diffraction on quantum chaos in non-superconducting systems was considered in ||. In NS 
structures the diffraction changes the orientation of kinetic momentum in a random fashion, 
which leads to spreading of the states over the whole system. We discuss manifesation of 
this effect below. 

For smooth disorder, one can derive the criterion of derealization via diffraction by using 
(H) and employing adiabatic approximation. For that, we make L in ([!]) position-dependent 
and interpret the energy ([[]) taken at p x = as spatially dependent potential energy. Also, 
we expand the square root in ([!]) and replace p x by —id/d x . This gives effective kinetic 
energy. Thus, one gets Hamiltonian: 

tt (n + |) v F f 1 5Ux) \ 

^effective - £ j ^ Q2 £ J W 

Here 5L(x) = L(x) — L is the deviation of the film thickness from its mean value. Since 
the particle is localized near the thickest place, where L(x) is maximal, one can write: 
5L(x) ~ A(l — x 2 /r1). Then one can estimate the width d of the ground state wave function 
by comparing (H) to the harmonic oscillator problem: 

Quantum effects do not destroy a localized state if its smearing given by ([T(]) is much less 
than the potential well width, i.e., if d r c . Thus, one arrives at another condition: 




(11) 

where r* is defined in (|8]). This is quantum criterion of localization. 

The classical and quantum criteria (§) and (|ll]) determine when localization can take 
place. Note that (||) and ( |TTD are compatible only when A > A^. This is expected, because 



if the thickness fluctuations are less than A^, the disorder cannot separate a group of states 
out of the continuum and localize these states. (In the case A < Xp, there is no room for 
an extra wavelength between the NS and metal-vacuum interfaces). 

Can one satisfy (|8p and flTT| ) in a real system? The film thickness L must be larger than, 
or of the order of the superconducting coherence length £ - Otherwise, due to the proximity 
gap induced in the normal layer there will be no excitations with the energies of interest. 
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Besides, the NS interface must have width > £o in order that the reflection is fully Andreev. 
Thus, for a superconductor with T c ~ 10K one gets L ~ lOOOA. Since Xp < A is required, 
let us take A ~ 10A. Thus, the criteria (|8]) and ( [TT]) give 10A < r c < lOOA. This means 
that the surface must have a certain degree of smoothness. At present, it is difficult to say 
how realistic this condition is. One can imagine a situation where all abrupt jumps of the 
surface are screened by conducting electrons so that the resulting potential is sufficiently 
smooth. 

Note that the localization of the type described above is quite different from the usual 
one. First, it occurs only if the scattering is classical, in contrast with the usual Aderson 
localization which is due to quantum nature of scattering. A manifestation of that is the 
suppression of localization in this system by short range disorder. Note, however, that at 
sufficiently high impurity concentration the electrons are again localized, now by Anderson 
mechanism. Thus, localization is reentrant with respect to disorder strength. We think that 
at high impurity concentration Andreev reflection should enhance Anderson localization 
effect. However, this question certainly needs more attention. 

Secondly, in this system there is no mobility edge: the energies of the localized and 
delocalized states are not separated. This apparently contradicts the standard ergodicity 
argument by Mott B about the absence of coexistence of localized and delocalized states 
with the same energy. The reason that there is no ergodicity [[H| in our problem is that the 
disorder is smooth. Due to this smoothness, there appear adiabatic barriers dividing the 
phase space into domains with very different dynamical characteristics (see Fig. |j). 

Finally, due to the presence of the superconductor, the localization in the NS struc- 
ture is less sensitive to Coulomb interaction effects. In fact, in this problem we deal with 
charged quasiparticles in a highly conducting and thus well-screening medium. Usually, the 
appearence of localized states at the metal-insulator transition in disordered systems is con- 
trolled by effects of Coulomb interaction. The reason for the importance of the interaction 
is that due to poor conductivity near the transition the screening of the interaction is very 
slow. In contrast, in the NS structure the interaction is screened by the superconductor. 
Therefore, the screening is always fast, no matter how slow individual localized electrons 
are. 

Finally, we discuss how localization of this unusual type could be observed. Perhaps, it 
cannot reveal itself through the conductivity, because the superconductor will always shunt 
electrical conductivity of the normal film. Instead, one can measure the local tunneling 
density of states, which could be probed, for example, by an STM. 

If the normal film is flat, the spectrum of electrons is described by ([I]). The average 
density of states corresponding to (|T|) has the well-known sawtooth structure. This structure 
should be independent of the STM tip position. 

If the localized states are present, they will add spatially dependent features to the 
tunneling density of states. Each localized state will give rise to a peak in the density of 
states if the STM tip is close to the place where the state is localized. Thus, one can study 
spatial correlation of peaks (and other features) in the local tunneling density of states. If 
the peaks are really due to localized states, they should be spatially uncorrelated. This 
measurement scheme is insensitive to the supercurrent. 

To summarize, we studied localization of quasiparticles in a normal metal film boarded 
by Andreev mirror. This system exhibits a new type of localization caused by self-retracing 
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due to Andreev reflection. We derived the criteria for this localization, and discussed its 
manifestation in the tunneling density of states. 
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FIG. 1. (a) Idealized picture of localization. The quasiparticle bounces between interfaces 
without spreading over the whole system, (b) At finite quasiparticle energy there is no perfect 
self-retracing in Andreev scattering. This results in a slow drift along the interface. 
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FIG. 2. Localization of the quasiparticle trajectory by surface fluctuation. The curvature of 
the surface eliminates average drift. 
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FIG. 4. Poincare section for y(x) = 10 + 1/(1 + x 2 ), e = W~ 3 ep. The central island and the 
four islands located symmetrically around it correspond to different periodic trajectories. If one 
denotes each collision with NS interface by S, and with metal-vacuum interface by V, then the 
central island corresponds to VSVSVS. . . sequence of collisions, while the lateral ones correspond 
to VVSVVS. . . sequence. 
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